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Abstract. We perform hydrodynamical simulations of core collapse supernovae (CCSNe) with 
. . . a cylindrically-symmetrical numerical code (FLASH) to study the inflation of bubbles and the 
initiation of the explosion within the frame of the jittering-jets model. We study the typical time- 
scale of the model and compare it to the typical time-scale of the delayed neutrino mechanism. 
Our analysis shows that the explosion energy of the delayed neutrino mechanism is an order of 
' ^ magnitude less than the required 10^^ erg. 
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1. The Jittering- Jets Model 

Accretion-outflow systems are commonly observed in astrophysics when compact ob- 
jects accrete mass via an accretion disk that launches jets. This is the basic engine 
(— I ' assumed in jet-based CCSN models (Papish & Soker 2011, Lazzati et al. 2011). Our 
Jittering- Jet Model for CCSN explosions is based on the following ingredients. (1) We do 
Q ' not try to revive the stalled shock. To the contrary, our model requires the material near 
the stalled-shock to fall inward and form an accretion disk around the newly born neu- 
■ tron star (NS) or black hole (BH). (2) We conjecture that due to stochastic processes and 
] the stationary accretion shock instability (SASI; e.g., Blondin & Shaw 2007, Fernandez 
2010) segments of the post-shock accreted gas possess local angular momentum (see also 
T— I ] Foglizzo et al. (2012) for an experimental demonstration). (3) We assume that the accre- 
^ ' tion disk launches two opposite jets. Due to the rapid change in the disk's axis, the jets 
can be intermittent and their direction rapidly varying. These are termed jittering jets. 
(4) The jets penetrate the infalling gas up to a distance of fewxlOOO km, i.e., beyond 
^vq the stalled-shock. The jets deposit their energy inside the star via shock waves, and form 
,• hot bubbles. (5) The jets are launched only in the last phase of accretion onto the NS. 

We perform 2.5D numerical simulations using the FLASH code (Fryxell et al. 2000). 
We use 2D cylindrical coordinates on a grid of size 1.5 x 10^ cm in each direction. We 
T-H use 10 mesh grid refinements which gives us a resolution of 3 km at the inner boundary 
^ of the grid of r = 70 km. For the initial conditions of the ambient gas in the core we 
•'^ \ used the ISA/© model of Liebendorfer et al. (2005) who made a ID simulations of the 
' core bounce. We start the simulation by launching jets at 0.25 s after bounce. We inject 
I a jet (only one jet as only one half of space is simulated) at 75 km from the center with 
" " " a full opening angle of 10 degrees. To simulate the jittering effect we inject the jet for a 
time interval of Atj — 0.05 s with a constant angle 6'„ relative to the z axis. After this 
time interval, we stop the jets for a period of Atp = 0.05 s. We then continue with a jet 
at a different angle of dn+i, for the same Atj = 0.05 s. We repeat this process for several 
times. The velocity of the jet is taken to be 10^° cm s~^ with with a mass outflow rate of 
M2j = 4 X 10^^ g s~^. The total energy carried by the two opposite jets combined over 
all episodes is i?2j = 2 x 10^^ erg. Our preliminary results are presented in Fig. 1. 
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Figure 1. Time t = 0.61 s after bounce witii tiie first jets episode starts at t — 0.25 s after 
bounce. 4 episodes of jets were launched up to tliis time, in directions 9 = 10° , 30° , 50° , and 
35°. The first two episodes are already mixed together; the 4th jet's episode was turned off 
0.01 seconds ago. Note that accretion from the equatorial plane (the horizontal plane in the 
figure) continues alongside with the outflow induced by the jets. The left plot shows the density, 
the right plot shows the log of the velocity. All values are in cgs units. The size of the grid is 
1.5 X 10^ cm in each direction. 

2. Time-scale Considerations 

The typical parameters of the model are the following. The active phase of the jets 
lasts for a typical time equals to few times the free fall time from the region where the 
baryonic mass is ^ 1.5 — 1.6Mq. This radius is ~ 3000 km and the free fall time is 0.4 s. 
The jets deposit their energy in a typical radius of rd — 10"* km. The average power of 
the jets is taken to be P2j — 10^^ erg s~^. The jets are active for ~ 1 — 2 s. This is 
also the dynamical time at 

U^^^ = l.6(-^f s, (2.1) 
Vcsc(rd) VlO-^km/ 

where v^sc is the escape velocity at r^, and where the central mass is taken to be I.SMq. 
The jets' power and the interaction time at ~ 10* km are consistent with the energy 
required to explode the star. From this time on, the hot bubbles (which might merge to 
one bubble) explode the star. 

We find that the same consideration of interaction time and explosion power is problem- 
atic for neutrino driven CCSN explosion models. In recent numerical results (e.g. Brandt 
et al. 2011, Hanke et al. 2011, Kuroda et al. 2012) the mechanical energy achieved by 
neutrino driven models is still significantly short of the desired ^ 10^^ erg required to 
explode the star. In the recent simulations of Mueller et al. (2012) the sum of the kinetic, 
thermal, and nuclear energy of the expanding gas in the core is a factor of ^ 4 smaller 
than the observed energy of CCSNe. We note that in Nordhaus et al. (2011) and Scheck 
et al. (2006) the explosion is achieved mainly by a continues wind. Here we refer to the 
delayed neutrino mechanism models where the energy of the wind is negligible. 

Let us apply the interaction time considerations to the delayed neutrino mechanism. 
The gain region, where neutrino heating is efficient, occurs in the region r ~ 200 — 700 km 
(Janka 2001). If energy becomes significant the gas will be accelerated and escape within a 
time of ^ r/[0.5tiosc (?")]■ From Janka (2001) we find the neutrino "optical depth" from r to 
infinity to be r ^ (r/100 km)~'^. The typical electron (and positron) neutrino luminosity 
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is Ljy ~ 5 X 10^^ erg s ^ (Mueller et al. 2012). Over all, if the interaction occurs near a 
radius r in the gain region, the energy that can be acquired by the expanding gas is 



We claim, therefore, that the total energy that can be used to revive the shock is limited 
to a typical value of ^ 10^° erg. This is along the recent results of numerical simulations 
of the delayed neutrino mechanism cited above. 

We conclude that the delayed neutrino explosion mechanism, where the explosion is 
due to neutrino heating in the gain region, as proposed by Wilson (Bethe & Wilson 1985), 
cannot work. It might lead to the reviving of the stalled shock under some circumstances, 
but it cannot lead to an explosion with an energy of 10^^ erg" 

Our basic conclusion is that if no ingredient is added to the neutrino delayed mech- 
anism, it falls short by an order of magnitude from the required energy to explode the 
star. Such an ingredient can be a strong wind, such as was applied by artificial energy 
deposition (Nordhaus et al. 2011). In their 2.5D simulations Scheck et al. (2006) achieved 
explosion that was mainly driven by a continues wind. The problem we see with winds is 
that they are less efficient than jets. Indeed, in order to obtain an explosion the winds in 
the simulations of Scheck et al. (2006) had to be massive. For that, in cases where they 
obtained energetic enough explosions the final mass of the NS was low (Mns < l.SM©). 
The problem we find is that such a wind must be active while accretion takes place; the 
accretion is required for the energy source. 

An infiow-outfiow situation naturally occurs with jets launched by accretion disks. For 
that, we propose the jittering-jet mechanism. Namely, we require the accretion process 
to continue for ~ 1 s. In our model there is no need to revive the accretion shock. 
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